Dynamic genome evolution in a model fern.

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

The biology of C. richardii as a tool to understand plant evolution.

The fern Ceratopteris richardii has been studied as a model organism for over 50 years because it is easy to grow and has a short life cycle. In particular, as the first homosporous vascular plant for which genomic resources were developed, C. richardii has been an important system for studying plant evolution. However, we know relatively little about the natural history of C. richardii. In this article, we summarize what is known about this aspect of C. richardii, and discuss how learning more about its natural history could greatly increase our understanding of the evolution of land plants.

Why Do Heterosporous Plants Have So Few Chromosomes?

The mechanisms controlling chromosome number, size, and shape, and the relationship of these traits to genome size, remain some of the least understood aspects of genome evolution. Across vascular plants, there is a striking disparity in chromosome number between homosporous and heterosporous lineages. Homosporous plants (comprising most ferns and some lycophytes) have high chromosome numbers compared to heterosporous lineages (some ferns and lycophytes and all seed plants). Many studies have investigated why homosporous plants have so many chromosomes. However, homospory is the ancestral condition from which heterospory has been derived several times. Following this phylogenetic perspective, a more appropriate question to ask is why heterosporous plants have so few chromosomes. Here, we review life history differences between heterosporous and homosporous plants, previous work on chromosome number and genome size in each lineage, known mechanisms of genome downsizing and chromosomal rearrangements, and conclude with future prospects for comparative research.

Cryptic diversity in the model fern genus Ceratopteris (Pteridaceae).

Cryptic species are present throughout the tree of life. They are especially prevalent in ferns, because of processes such hybridization, polyploidy, and reticulate evolution. In addition, the simple morphology of ferns limits phenotypic variation and makes it difficult to detect cryptic species. The model fern genus Ceratopteris has long been suspected to harbor cryptic diversity, in particular within the highly polymorphic C. thalictroides. Yet no studies have included samples from throughout its pan-tropical range or utilized genomic sequencing, making it difficult to assess the full extent of cryptic variation within this genus. Here, we present the first multilocus phylogeny of the genus using reduced representation genomic sequencing (RADseq) and examine population structure, phylogenetic relationships, and ploidy level variation. We recover similar species relationships found in previous studies, find support for the cryptic species C. gaudichaudii as genetically distinct, and identify novel genomic variation within two of the mostly broadly distributed species in the genus, C. thalictroides and C. cornuta. Finally, we detail the utility of our approach for working on cryptic, reticulate groups of ferns. Specifically, it does not require a reference genome, of which there are very few available for ferns. RADseq is a cost-effective way to work with study groups lacking genomic resources, and to obtain the thousands of nuclear markers needed to untangle species complexes.

How Many Tree Species of Birch Are in Alaska? Implications for Wetland Designations.

Wetland areas are critical habitats, especially in northern regions of North America. Wetland classifications are based on several factors, including the presence of certain plant species and assemblages of species, of which trees play a significant role. Here we examined wetland species of birch (Betula) in North America, with a focus on Alaska, and the use of birche tree species in wetland delineation. We sampled over 200 trees from sites, including Alaska, Alberta, Minnesota, and New Hampshire. We used genetic data from over 3000 loci detected by restriction site associated DNA analysis. We used an indirect estimate of ploidy based on allelic ratios and we also examined population genetic structure. We find that inferred ploidy is strongly associated with genetic groupings. We find two main distinct groups; one found throughout most of Alaska, extending into Alberta. This group is probably attributable to Betula kenaica, Betula neoalaskana, or both. This group has a diploid genetic pattern although this could easily be a function of allopolyploidy. The second major genetic group appears to extend from Eastern North America into parts of southeastern Alaska. This group represents Betula papyrifera, and is not diploid based on allelic ratios. Published chromosome counts indicate pentaploidy. Because B. papyrifera is the only one of the above species that is distinctly associated with wetland habitats, our findings indicate that tree species of birch found in most parts of Alaska are not reliable indicators of wetland habitats. These results help to support stronger wetland ratings assigned to the tree species of birch for delineation purposes.

The C-Fern (Ceratopteris richardii) genome: insights into plant genome evolution with the first partial homosporous fern genome assembly.

Ferns are notorious for possessing large genomes and numerous chromosomes. Despite decades of speculation, the processes underlying the expansive genomes of ferns are unclear, largely due to the absence of a sequenced homosporous fern genome. The lack of this crucial resource has not only hindered investigations of evolutionary processes responsible for the unusual genome characteristics of homosporous ferns, but also impeded synthesis of genome evolution across land plants. Here, we used the model fern species Ceratopteris richardii to address the processes (e.g., polyploidy, spread of repeat elements) by which the large genomes and high chromosome numbers typical of homosporous ferns may have evolved and have been maintained. We directly compared repeat compositions in species spanning the green plant tree of life and a diversity of genome sizes, as well as both short- and long-read-based assemblies of Ceratopteris. We found evidence consistent with a single ancient polyploidy event in the evolutionary history of Ceratopteris based on both genomic and cytogenetic data, and on repeat proportions similar to those found in large flowering plant genomes. This study provides a major stepping-stone in the understanding of land plant evolutionary genomics by providing the first homosporous fern reference genome, as well as insights into the processes underlying the formation of these massive genomes.

Worldwide relationships in the fern genus Pteridium (bracken) based on nuclear genome markers.

PREMISE: Spore-bearing plants are capable of dispersing very long distances. However, it is not known if gene flow can prevent genetic divergence in widely distributed taxa. Here we address this issue, and examine systematic relationships at a global geographic scale for the fern genus Pteridium. METHODS: We sampled plants from 100 localities worldwide, and generated nucleotide data from four nuclear genes and two plastid regions. We also examined 2801 single nucleotide polymorphisms detected by a restriction site-associated DNA approach. RESULTS: We found evidence for two distinct diploid species and two allotetraploids between them. The "northern" species (Pteridium aquilinum) has distinct groups at the continental scale (Europe, Asia, Africa, and North America). The northern European subspecies pinetorum appears to involve admixture among all of these. A sample from the Hawaiian Islands contained elements of both North American and Asian P. aquilinum. The "southern" species, P. esculentum, shows little genetic differentiation between South American and Australian samples. Components of African genotypes are detected on all continents. CONCLUSIONS: We find evidence of distinct continental-scale genetic differentiation in Pteridium. However, on top of this is a clear signal of recent hybridization. Thus, spore-bearing plants are clearly capable of extensive long-distance gene flow; yet appear to have differentiated genetically at the continental scale. Either gene flow in the past was at a reduced level, or vicariance is possible even in the face of long-distance gene flow.

ReFernment: An R package for annotating RNA editing in plastid genomes.

PREMISE OF THE STUDY: In the absence of cDNA, the annotation of RNA editing in plastomes must be done manually, representing a significant time cost to those studying the organellar genomes of ferns and hornworts. METHODS AND RESULTS: We developed an R package to automatically annotate apparent nonsense mutations in plastid genomes. The software successfully annotates such sites and results in no false positives for data with no sequencing or assembly errors. CONCLUSIONS: Compared to manual annotation, ReFernment offers greater speed and accuracy for annotating RNA editing sites. This software should be especially useful for researchers generating large numbers of plastome sequences for taxa with high levels of RNA editing.

Mobile Elements Shape Plastome Evolution in Ferns.

Plastid genomes display remarkable organizational stability over evolutionary time. From green algae to angiosperms, most plastid genomes are largely collinear, with only a few cases of inversion, gene loss, or, in extremely rare cases, gene addition. These plastome insertions are mostly clade-specific and are typically of nuclear or mitochondrial origin. Here, we expand on these findings and present the first family-level survey of plastome evolution in ferns, revealing a novel suite of dynamic mobile elements. Comparative plastome analyses of the Pteridaceae expose several mobile open reading frames that vary in sequence length, insertion site, and configuration among sampled taxa. Even between close relatives, the presence and location of these elements is widely variable when viewed in a phylogenetic context. We characterize these elements and refer to them collectively as Mobile Open Reading Frames in Fern Organelles (MORFFO). We further note that the presence of MORFFO is not restricted to Pteridaceae, but is found across ferns and other plant clades. MORFFO elements are regularly associated with inversions, intergenic expansions, and changes to the inverted repeats. They likewise appear to be present in mitochondrial and nuclear genomes of ferns, indicating that they can move between genomic compartments with relative ease. The origins and functions of these mobile elements are unknown, but MORFFO appears to be a major driver of structural genome evolution in the plastomes of ferns, and possibly other groups of plants.

Target sequence capture of nuclear-encoded genes for phylogenetic analysis in ferns.

PREMISE OF THE STUDY: Until recently, most phylogenetic studies of ferns were based on chloroplast genes. Evolutionary inferences based on these data can be incomplete because the characters are from a single linkage group and are uniparentally inherited. These limitations are particularly acute in studies of hybridization, which is prevalent in ferns; fern hybrids are common and ferns are able to hybridize across highly diverged lineages, up to 60 million years since divergence in one documented case. However, it not yet clear what effect such hybridization has on fern evolution, in part due to a paucity of available biparentally inherited (nuclear-encoded) markers. METHODS: We designed oligonucleotide baits to capture 25 targeted, low-copy nuclear markers from a sample of 24 species spanning extant fern diversity. RESULTS: Most loci were successfully sequenced from most accessions. Although the baits were designed from exon (transcript) data, we successfully captured intron sequences that should be useful for more focused phylogenetic studies. We present phylogenetic analyses of the new target sequence capture data and integrate these into a previous transcript-based data set. DISCUSSION: We make our bait sequences available to the community as a resource for further studies of fern phylogeny.

Fern genomes elucidate land plant evolution and cyanobacterial symbioses.

Ferns are the closest sister group to all seed plants, yet little is known about their genomes other than that they are generally colossal. Here, we report on the genomes of Azolla filiculoides and Salvinia cucullata (Salviniales) and present evidence for episodic whole-genome duplication in ferns-one at the base of 'core leptosporangiates' and one specific to Azolla. One fern-specific gene that we identified, recently shown to confer high insect resistance, seems to have been derived from bacteria through horizontal gene transfer. Azolla coexists in a unique symbiosis with N2-fixing cyanobacteria, and we demonstrate a clear pattern of cospeciation between the two partners. Furthermore, the Azolla genome lacks genes that are common to arbuscular mycorrhizal and root nodule symbioses, and we identify several putative transporter genes specific to Azolla-cyanobacterial symbiosis. These genomic resources will help in exploring the biotechnological potential of Azolla and address fundamental questions in the evolution of plant life.

Admixture, evolution, and variation in reproductive isolation in the Boechera puberula clade.

BACKGROUND: Hybridization is very common in plants, and the incorporation of new alleles into existing lineages (i.e. admixture) can blur species boundaries. However, admixture also has the potential to increase standing genetic variation. With new sequencing methods, we can now study admixture and reproductive isolation at a much finer scale than in the past. The genus Boechera is an extraordinary example of admixture, with over 400 hybrid derivates of varying ploidy levels. Yet, few studies have assessed admixture in this genus on a genomic scale. RESULTS: In this study, we used Genotyping-by-Sequencing (GBS) to clarify the evolution of the Boechera puberula clade, whose six members are scattered across the western United States. We further assessed patterns of admixture and reproductive isolation within the group, including two additional species (B. stricta and B. retrofracta) that are widespread across North America. Based on 14,815 common genetic variants, we found evidence for some cases of hybridization. We find evidence of both recent and more ancient admixture, and that levels of admixture vary across species. CONCLUSIONS: We present evidence for a monophyletic origin of the B. puberula group, and a split of B. puberula into two subspecies. Further, when inferring reproductive isolation on the basis of presence and absence of admixture, we found that the accumulation of reproductive isolation between species does not seem to occur linearly with time since divergence in this system. We discuss our results in the context of sexuality and asexuality in Boechera.

Serologic Evidence of Various Arboviruses Detected in White-Tailed Deer (Odocoileus virginianus) in the United States.

White-tailed deer (Odocoileus virginianus) are an abundant mammal with a wide geographic distribution in the United States, which make them good sentinels for monitoring arboviral activity across the country. Exposure to various arboviruses has been detected in white-tailed deer, typically in conjunction with another diagnostic finding. To better assess the exposure of white-tailed deer to seven arboviruses, we tested 1,508 sera collected from 2010 to 2016 for antibodies to eastern equine encephalitis (2.5%), Powassan (4.2%), St. Louis encephalitis, (3.7%), West Nile (6.0%), Maguari (19.4%), La Crosse (30.3%), and bluetongue (7.8%) viruses. At least one arbovirus was detected in 51.3%, and exposure to more than one arbovirus was identified in 17.6% of the white-tailed deer sampled.

A genetic signature of the evolution of loss of flight in the Galapagos cormorant.

We have a limited understanding of the genetic and molecular basis of evolutionary changes in the size and proportion of limbs. We studied wing and pectoral skeleton reduction leading to flightlessness in the Galapagos cormorant (Phalacrocorax harrisi). We sequenced and de novo assembled the genomes of four cormorant species and applied a predictive and comparative genomics approach to find candidate variants that may have contributed to the evolution of flightlessness. These analyses and cross-species experiments in Caenorhabditis elegans and in chondrogenic cell lines implicated variants in genes necessary for transcriptional regulation and function of the primary cilium. Cilia are essential for Hedgehog signaling, and humans affected by skeletal ciliopathies suffer from premature bone growth arrest, mirroring skeletal features associated with loss of flight.

Surveillance for Highly Pathogenic Avian Influenza in Wild Turkeys ( Meleagris gallopavo ) of Minnesota, USA during 2015 Outbreaks in Domestic Poultry.

An outbreak of a novel reassortant of highly pathogenic avian influenza A (H5N2) virus (HPAIV) decimated domestic turkeys ( Meleagris gallopavo ) from March through mid-June, 2015 in the state of Minnesota, US. In response, as part of broader surveillance efforts in wild birds, we designed a pilot effort to sample and test hunter-harvested Wild Turkeys ( Meleagris gallopavo ) for HPAIV in Minnesota counties with known infected poultry facilities. We also collected opportunistic samples from dead Wild Turkeys or live Wild Turkeys showing neurologic signs (morbidity and mortality samples) reported by the public or state agency personnel. Cloacal and tracheal samples were collected from each bird and screened for avian influenza virus (AIV) RNA by real-time reverse transcription PCR. From 15 April to 28 May 2015, we sampled 84 hunter-harvested male Wild Turkeys in 11 Minnesota counties. From 7 April 2015 through 11 April 2016, we sampled an additional 23 Wild Turkeys in 17 Minnesota counties. We did not detect type A influenza or HPAIV from any samples, and concluded, at the 95% confidence level, that apparent shedding prevalence in male Wild Turkeys in central Minnesota was between 0% and 2.9% over the sampling period. The susceptibility of wild turkeys to HPAIV is unclear, but regular harvest seasons make this wild gallinaceous bird readily available for future AIV testing.

Evolutionary Conservation of ABA Signaling for Stomatal Closure.

Abscisic acid (ABA)-driven stomatal regulation reportedly evolved after the divergence of ferns, during the early evolution of seed plants approximately 360 million years ago. This hypothesis is based on the observation that the stomata of certain fern species are unresponsive to ABA, but exhibit passive hydraulic control. However, ABA-induced stomatal closure was detected in some mosses and lycophytes. Here, we observed that a number of ABA signaling and membrane transporter protein families diversified over the evolutionary history of land plants. The aquatic ferns Azolla filiculoides and Salvinia cucullata have representatives of 23 families of proteins orthologous to those of Arabidopsis (Arabidopsis thaliana) and all other land plant species studied. Phylogenetic analysis of the key ABA signaling proteins indicates an evolutionarily conserved stomatal response to ABA. Moreover, comparative transcriptomic analysis has identified a suite of ABA-responsive genes that differentially expressed in a terrestrial fern species, Polystichum proliferum These genes encode proteins associated with ABA biosynthesis, transport, reception, transcription, signaling, and ion and sugar transport, which fit the general ABA signaling pathway constructed from Arabidopsis and Hordeum vulgare The retention of these key ABA-responsive genes could have had a profound effect on the adaptation of ferns to dry conditions. Furthermore, stomatal assays have shown the primary evidence for ABA-induced closure of stomata in two terrestrial fern species Pproliferum and Nephrolepis exaltata In summary, we report, to our knowledge, new molecular and physiological evidence for the presence of active stomatal control in ferns.

Surveillance for highly pathogenic H5 avian influenza virus in synanthropic wildlife associated with poultry farms during an acute outbreak.

In November 2014, a Eurasian strain H5N8 highly pathogenic avian influenza virus was detected in poultry in Canada. Introduced viruses were soon detected in the United States and within six months had spread to 21 states with more than 48 million poultry affected. In an effort to study potential mechanisms of spread of the Eurasian H5 virus, the United States Department of Agriculture coordinated several epidemiologic investigations at poultry farms. As part of those efforts, we sampled synanthropic birds and mammals at five infected and five uninfected poultry farms in northwest Iowa for exposure to avian influenza viruses. Across all farms, we collected 2,627 samples from 648 individual birds and mammals. House mice were the most common mammal species captured while house sparrows, European starlings, rock pigeons, swallows, and American robins were the most commonly captured birds. A single European starling was positive for Eurasian H5 viral RNA and seropositive for antibodies reactive to the Eurasian H5 virus. Two American robins were also seropositive. No mammal species showed evidence of infection. These results indicate synanthropic species merit further scrutiny to better understand potential biosecurity risks. We propose a set of management practices aimed at reducing wildlife incursions.

Surveillance for Highly Pathogenic Avian Influenza Virus in Wild Birds during Outbreaks in Domestic Poultry, Minnesota, 2015.

In 2015, a major outbreak of highly pathogenic avian influenza virus (HPAIV) infection devastated poultry facilities in Minnesota, USA. To understand the potential role of wild birds, we tested 3,139 waterfowl fecal samples and 104 sick and dead birds during March 9-June 4, 2015. HPAIV was isolated from a Cooper's hawk but not from waterfowl fecal samples.

Antibody Detection and Molecular Characterization of Toxoplasma gondii from Bobcats (Lynx rufus), Domestic Cats (Felis catus), and Wildlife from Minnesota, USA.

Little is known of the epidemiology of toxoplasmosis in Minnesota. Here, we evaluated Toxoplasma gondii infection in 50 wild bobcats (Lynx rufus) and 75 other animals on/near 10 cattle farms. Antibodies to T. gondii were assayed in serum samples or tissue fluids by the modified agglutination test (MAT, cut-off 1:25). Twenty nine of 50 bobcats and 15 of 41 wildlife trapped on the vicinity of 10 farms and nine of 16 adult domestic cats (Felis catus) and six of 14 domestic dogs resident on farms were seropositive. Toxoplasma gondii oocysts were not found in feces of any felid. Tissues of all seropositive wild animals trapped on the farm were bioassayed in mice and viable T. gondii was isolated from two badgers (Taxidea taxus), two raccoons (Procyon lotor), one coyote (Canis latrans), and one opossum (Didelphis virginiana). All six T. gondii isolates were further propagated in cell culture. Multi-locus PCR-RFLP genotyping using 10 markers (SAG1, SAG2 (5'-3'SAG2, and alt.SAG2), SAG3, BTUB, GRA6, c22-8, c29-2, L358, PK1, and Apico), and DNA from cell culture derived tachyzoites revealed three genotypes; #5 ToxoDataBase (1 coyote, 1 raccoon), #1 (1 badger, 1 raccoon, 1 opossum), and #2 (1 badger). This is the first report of T. gondii prevalence in domestic cats and in bobcats from Minnesota, and the first isolation of viable T. gondii from badger.

Origins and diversity of a cosmopolitan fern genus on an island archipelago.

Isolated oceanic islands are characterized by patterns of biological diversity different from that on nearby continental mainlands. Isolation can provide the opportunity for evolutionary divergence, but also set the stage for hybridization between related taxa arriving from different sources. Ferns disperse by haploid spores, which are produced in large numbers and can travel long distances in air currents, enabling these plants to become established on most oceanic islands. Here, we examine the origins and patterns of diversity of the cosmopolitan fern genus Pteridium (Dennstaedtiaceae; bracken) on the Galapagos Islands. We use nucleotide sequences from two plastid genes, and two nuclear gene markers, to examine phylogeography of Pteridium on the Galapagos Islands. We incorporate data from a previous study to provide a worldwide context. We also sampled new specimens from South and Central America. We used flow cytometry to estimate genome size of some accessions. We found that both plastid and nuclear haplotypes fall into two distinct clades, consistent with a two-diploid-species taxonomy of P. aquilinum and P. esculentum. As predicted, the allotetraploid P. caudatum possesses nuclear haplotypes from both diploid species. Samples from the Galapagos include P. esculentum subsp. arachnoideum, P. caudatum and possible hybrids between them. Multiple Pteridium taxa were also observed growing together at some sites. We find evidence for multiple origins of Pteridium on the Galapagos Islands and multiple origins of tetraploid P. caudatum throughout its range in Central and South America. We also posit that P. caudatum may include recent diploid hybrids, backcrosses to P. esculentum, as well as allotetraploid plants. The Galapagos Islands are positioned close to the equator where they can receive dispersing propagules from both hemispheres. This may partly explain the high levels of diversity found for this cosmopolitan fern on these islands.